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ABSTRACT 

In the ongoing HATNet survey we have detected a giant planet, with radius 1.33±0.06 i?jup and mass 
1.06 ± 0.12 Mjup, transiting the bright {V = 10.5) star GSC 03239-00992. The planet is in a circular 
orbit with period 3.852985 ± 0.000005 days and mid-transit epoch 2,454,035.67575 ± 0.00028 (HJD). 
The parent star is a late F star with mass 1.29 ± 0.06 M©, radius 1.46 ± 0.06 Rq, T^s ^ 6570 ± 80 K, 
[Fe/H] = -0.13 ± 0.08 and age ~ 2.3lo;7Gy. With this radius and mass, HAT-P-6b has somewhat 
larger radius than theoretically expected. We describe the observations and their analysis to determine 
physical properties of the HAT-P-6 system, and briefly discuss some implications of this finding. 
Subject headings: stars: individual (GSC 03239-00992, HAT-P-6) - planetary systems 



1. INTRODUCTION 

The detection of transiting exoplanets is very impor- 
tant to exoplanet research because of the information 
about both planetary radius and mass that comes from 
photometric transit light curves combined with follow- 
up radial velocity observations. The transiting exoplan- 

' ets known as of this writing span a wide range in the 
physical parameter space of planetary mass, radius, or- 
bital period, semi-major axis, eccentricity; and parent 

I star parameters, including mass, radius, effective tem- 
perature, metallicity, and age. Filling out their distri- 
bution in this multidimensional space is certain to give 
us important information on the origin and evolution of 
exoplanetary systems. Here we report on the discovery 

■ by HATNet of its sixth transiting planet, HAT-P-6b, an 
inflated Jupiter-mass gas giant in an essentially circu- 
lar orbit about an F dwarf star with slightly sub-solar 
metallicity. 

2. PHOTOMETRIC DETECTION 

! The HATNet telescopes HAT-6 and HAT-9 (HATNet; 
' Bakos et al. 2002, 2004) observed HATNet field G161, 
, centered at a = 00'^32"\ S = -|-37°30', on a near-nightly 
basis from 2005 August 12 to 2005 December 16. We 
have gathered altogether 9550 5-min exposures, each 
yielding photometric measurements for approximately 
35,000 stars with / < 14 and about 10,000 stars with 
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better than 2% light curve rms. The field was observed 
in network mode, whereby at the end of its nightly ob- 
serving sequence the HAT-6 telescope in Arizona handed 
off to the HAT-9 telescope in Hawaii, thus extending the 
duration of continuous observations. Following standard 
frame calibration procedures, astrometry was performed 
as described in Pal & Bakos (2006), and aperture pho- 
tometry results were subjected to External Parameter 
Decorrelation (EPD, described briefly in Bakos et al. 
2007a), and the Trend Filtering Algorithm (TEA; Kovacs 
et al. 2005). We searched the light curves for box-shaped 
transit signals using the BLS algorithm of Kovacs et al. 
(2002). A very significant periodic dip in intensity was 
detected in the / « 10.6 magnitude star GSC 03239- 
00992 (also known as 2MASS 23390581+4227575, with 
a depth of 9 mmag, a period of P = 3.8529 days, and a 
duration of 3.1 hours. 
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Fig. 1. — (a) Unbinned instrumental 7-band discovery light 
curve of HAT-P-6 obtained with HATNet, folded with the pe- 
riod of P = 3.852985 days. Superimposed (larger dots) is the 
same data binned to 1/200 in phase, (b) Unbinned instrumen- 
tal Sloan 2-band photometry collected with the KeplerCam at the 
FLWO 1.2 m telescope on 2006 October 26 (top curve) and again 
on UT 2007 September 4 (next curve); superimposed on both is 
our best- fit transit model curve (see text). 

The HATNet discovery light curve is shown in Fig. la. 
As is shown in the following sections we deduce that the 
signal is due to the transit of a Jovian planet across the 
face of the star. Hereafter we refer to the star as HAT- 
P-6, and to the planetary companion as HAT-P-6b. 
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TABLE 1 

Relative radial velocity measurements of HAT-P-6 
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3. FOLLOW-UP OBSERVATIONS 

HAT-P-6 was observed spectroscopically with the CfA 
Digital Speedometer (Latham 1992) at the FLWO 1.5 m 
Tilhnghast reflector of the Fred L. Whipple Observatory 
(FLWO) in order to rule out the possibility that the ob- 
served drop in brightness is caused by a transiting low- 
mass stellar companion rather than a planet, as well as to 
characterize the rotation and surface gravity of the star. 
Seven spectra were obtained over an interval of 92 days. 
Radial velocities were obtained by cross-correlation and 
have a typical precision of 0.4 km s^^. They showed no 
variation within the uncertainties, ruling out a compan- 
ion of stellar mass. The mean heliocentric radial velocity 
is -22.7± O.Skms-^ 

Photometric follow-up of HAT-P-6 was then carried 
out in the Sloan 2;-band with KeplerCam (see e.g. Hol- 
man et al. 2007) on the FLWO 1.2 m telescope, on 2006 
October 26. An astrometric solution between the indi- 
vidual frames and the 2MASS catalog was carried out 
using first order polynomials based on ~400 stars per 
frame. Aperture photometry was performed using a se- 
ries of apertures in fixed positions around the 2MASS- 
bascd {x, y) pixel coordinates. We selected a frame taken 
near the meridian and used ~260 stars and their magni- 
tudes as measured on this reference frame to transform 
all other frames to a common instrumental magnitude 
system. The aperture yielding the lowest scatter outside 
of transit was used in the subsequent analysis. The light 
curve was then de-correlated against trends using the 
out-of-transit sections and a dependence on hour angle 
(see Fig. lb, upper curve). 

A follow-up KeplerCam observation was recently ob- 
tained (2007 Sep 4), for the purpose of improving the 
photometric accuracy of the transit curve model, and 
also to determine the orbital period and mid- transit time 
with maximum accuracy. The data were treated identi- 
cally to those for the 2006 October 26 transit; results are 
shown in the lower curve of Fig. lb. 

Following the first KeplerCam observation, high res- 
olution spectroscopy was initiated with the HIRES in- 
strument (Vogt et al. 1994) on the Keck I telescope, in 
order both to determine the stellar parameters more pre- 
cisely and to characterize the radial velocity signal due 
to the companion. With a spectrometer slit of 0'.'86 the 
resolving power is A/AA w 55,000, and the wavelength 



coverage is ^ 3800 — 8000 A. An iodine gas absorption 
cell was used to superimpose a dense forest of I2 lines 
on the stellar spectrum and establish a highly accurate 
wavelength fiducial (see Marcy & Butler 1992). In to- 
tal 15 exposures were obtained between 2006 October 14 
and 2007 August 25 with the iodine cell, along with one 
without I2 for use as a template. Relative radial veloci- 
ties in the Solar System barycentric frame were derived 
as described by Butler et al. (1996), incorporating full 
modeling of the spatial and temporal variations of the 
instrumental profile. Data and their internal errors are 
listed in Table 1. 

4. ANALYSIS 

We determined the parameters of the star, and of the 
transiting planet, from the combined photometric and 
spectroscopic data by the following procedure. 

First, the iodine-free template spectrum from Keck was 
used for an initial determination of the atmospheric pa- 
rameters of the star. Spectral synthesis modeling was 
carried out using the SME software (Valenti & Piskunov 
1996), with wavelength ranges and atomic line data as 
described by Valenti & Fischer (2005). We obtained 
initial values as follows: effective temperature T^e = 
6353 ± 88K, surface gravity logg* = 3.84 ± 0.12, iron 
abundance [Fe/H] = —0.23 ± 0.08, and projected rota- 
tional velocity wsin « = 8.7 ± 1.0 kms^^. The tempera- 
ture and surface gravity correspond to an F8 dwarf. The 
uncertainties in SME-derived parameters quoted here 
and in the remainder of this discussion are twice the sta- 
tistical uncertainties; this reflects our attempt, based on 
prior experience, to incorporate systematic errors. 

Next, an initial modeling of the 2006 KeplerCam light 
curve was carried out using the formalism based on 
Mandel & Agol (2002), using the quadratic model for 
limb darkening without the assumption that the planet 
is small. The initial quadratic limb-darkening coefficients 
(namely, 71 = 0.1452 and 72 — 0.3488) were taken from 
the tables of Claret (2004) by interpolation to the above- 
mentioned SME values. The period was initially fixed at 
the value given earlier from the HATNet photometry, and 
the orbit was assumed to be circular, based on the results 
from the radial velocity fit to be discussed below. The 
four adjustable parameters are the planet-to-star ratio 
of the radii (i?p/i?*), the normalized separation (a/i?*) 
where a is the semi-major axis of the relative orbit, the 
normalized impact parameter (6 = acosi/i?*), and the 
time of the center of the transit (Tc). Both for the ini- 
tial light curve fit and the further fits (sec below) we used 
the Markov Chain Monte-Carlo (MCMC) method to find 
the best fit parameters (see, e.g. Ford 2004, for a com- 
prehensive description). To estimate errors we used the 
method of refitting to synthetic data sets to determine 
their uncertainties and correlations. The synthetic data 
sets consist of the analytic (fitted) model plus Gaussian 
and red-noise, based on both white-noise and red-noise 
estimates of the residuals. The characteristics of the red 
noise component of the residuals were preserved by per- 
turbing only the phase of their Fourier spectrum. We 
found that this method of error estimation is preferable 
to direct estimation using the MCMC method, since it 
is not sensitive to the number of out-of-transit points 
used. We also note that instead of a/ii* and 6, we used 
the parameters C/R* and b"^ for the fit, where ( is an 
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TABLE 2 
Stellar parameters for HAT-P-6. 



Parameter Value Source 



Tgff (K) 6570 ± 80 SME'' 

[Fe/H] -0.13 ± 0.08 SME 

logg* (cgs) 4.22 ±0.03 Y^+LC+SME^ 

i;sini (kms-i) .... 8.7±1.0 SME 

M* (Mg) 1.29 ±0.06 Y2+LC+SME 

(Rq) 1.46 ±0.06 Y2+LC+SME 

(Lq) 3.57l°-g Y^+LC+SME 

Mv (mag) 3.36 ±0.16 Y^+LC+SME 

Age (Gyr) 2.31^1 Y^+LC+SME 

Distance (pc) 260 ± 20 Y^+LC+SME 



'^SME = Package for analysis of high-resolution spectra Valenti 
& Piskunov (1996). 

bY2+LC+SME = Yale-Yonsei isochrones (Yi et al. 2001), light 
curve parameters, and SME results. 

auxiliary variable with dimensions of velocity, defined by 
( = 27ra/ (PVI — b'^) ■ These parameters have been cho- 
sen to eliminate the correlation between a/i?* and 6^ (see 
also Bakos et al. 2007b). 

Next, we used the values of T^g and [Fe/H] from the 
initial SME analysis, together with the initial value of 
a/Ri, from the Mandcl-Agol fit, to estimate the stellar 
properties from comparison with the Yonsei-Yale (Y^) 
stellar evolution models by Yi et al. (2001). The value of 
a/R^, is closely related to the stellar density, and is thus a 
proxy for luminosity (i*). As described by Sozzetti et al. 
(2007), a/Ri, is typically a better constraint on L^, than 
the spectroscopic value of logg*, which has a relatively 
subtle effect on the line profiles and whose determina- 
tion is therefore more susceptible to systematic errors. 
Following Sozzetti et al. (2007) we determined the range 
of stellar masses and radii that arc consistent with the 
SME-determined values ofTes, [Fe/H], and a/i?* derived 
from the light curve. We obtained = l.lQtg l^Mo 
and Ri, = 1A5'^q1\ Rq. The resultant surface gravity, 

log 5* = 4.19to.ioi "^^s significantly larger than the ini- 
tial SME-derived value discussed above. We fixed log g^ 
at this new, more accurate value, and repeated the SME 
analysis. Fixing surface gravity slightly affects the so- 
lution for other parameters such as [Fe/H] and other 
individual element abundances as well as vsini. The 
correlations between these free parameters in turn affect 
the final solution for effective temperature, causing it to 
change by about twice its original uncertainty. The re- 
sulting values from this iteration are Teff = 6570 ± 80 K, 
[Fe/H] = -0.13±0.08, and vsini = 8.2 ± LOkms"!. 

We then performed a second iteration of the above 
steps using the new values of Toff, log g*, and [Fe/H], as 
well as correspondingly changed limb darkening coeffi- 
cents: 71 = 0.1211 and 72 = 0.3646. This iteration also 
incorporated both the 2006 and 2007 KeplerCam light 
curves (Fig. lb) into a single fit, which yielded the times 
of each transit center Td and Tc2, as well as the shape 
parameters a/R^,, Rp/R^,, and 6, which we assumed to 
be the same for both of the transit events. 

Finally, using the new values of a/i?*, Teff, and [Fe/H] 
together with the Y^ evolutionary models discussed 
above, we re-determined the stellar parameters. Table 2 
summarizes these, which include a further refined deter- 
mination of log 51* = 4.217lo!o27- This value is the same 



TABLE 3 

Spectroscopic and light curve solutions for HAT-P-6, and 

INFERRED PLANET PARAMETERS 



Parameter 


Value 


Light curve parameters 




P (days) 


3.852985 ± 0.000005 


Tc (HJD-2,400,000) 


54,035.67575 ± 0.00028 


Tm (days)^ 


0.1461 ± 0.0017 


-tl2 — -^34 t^cidybj 


U.U-LOO in U.UUl-L 


a/R^ 


7.69 ± 0.22 




0.09338 ± 0.00053 


b = a, cos i j Ti,^ 


0.602 ± 0.030 


% (deg) 


85? 51 ± 0?35 


C (i?*/day) 


15.696 ± 0.086 


Spectroscopic parameters 




K (ms-i) 


115.5 ±4.2 


7 (kms~^) 


-22.7 ±0.5 


e 


(adopted) 


Planet parameters 




(Mj) 


1.057 ±0.119 


Rv (R}) 


1.330 ± 0.061 


Pp (gcm-3) 


0.558 ± 0.047 


a (AU) 


0.05235 ± 0.00087 


Qp (ms-2) 


14.8 ± 1.0 



* T14: total transit duration, time between first to last contact; 
T12 = T34: ingress/egress time, time between first and second, or 
third and fourth contact. 



as the one in the previous iteration to within uncertain- 
ties, so we stopped the iteration at this point. 

^From the value of vsini and the stellar radius ii*, 
and assuming an inclination of the rotational equator of 
approximately 90 deg (that is, similar to the inclination 
of the planetary orbital plane), we estimate the stellar 
rotation period to be about 9 days. We have searched 
the HATNet light curve (after elimination of data during 
transits) for periodicities between 6 and 12 days, such as 
might be produced by significant spots or other stellar 
activity, and found no evidence for this. 

An independent estimate of the age was obtained from 
the Ca"*" H and K line emission strength, measured from 
the mean of 14 Keck spectra: log Rhk = —4.81, lead- 
ing to an age of 2.8 Gy (based on Noyes et al. 1984), 
consistent with the evolutionary track age. The lu- 
minosity and Teff, imply an absolute visual magnitude 
My = 3.36 ± 0.16 mag. Combined with the apparent vi- 
sual magnitude {V = 10.440 ± 0.036; Droege et al. 2006; 
Hog et al. 2000), this yields a distance of 260 ± 20 pc, ig- 
noring extinction. For reference, the final values for the 
stellar parameters listed above are given in Table 2. 

Setting the difference between the transit centers of 
the two KeplerCam light curves, Tc2 — Td, to 81 or- 
bital periods, we obtain a more accurate value for the 
period: P = 3.852985 ± 0.000005 days. This value, 
and the reference epoch of mid-transit, Tc = Td = 
2,454,035.67575 ± 0.00028, are given in Table 3. Also 
given in Table 3 are the final fitted light curve parameters 
a/Ri,, Rp/Ri,, b, i (orbital inclination), and the radius of 
the planet, Rp = 1.330 ± 0.061 as determined from 
Ri, and Rp/R^,. The modeled z-band transit light curve 
is shown as a solid line superimposed on the data points 
in Fig. lb. 

The Keck radial velocity data were initially fit with a 
Keplerian model with period and epoch constrained to 
the value determined from the photometry (§2) and no 
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Fig. 2. — (a) Radial- velocity measurements from Keck for HAT- 
P-6, along with an orbital fit, shown as a function of orbital phase. 
The center-of-mass velocity has been subtracted, (b) Phased resid- 
uals after subtracting off the orbital fit. The rms variation of the 
residuals is about 8.6m s^-'. (c) Bisector spans (BS) for the 16 
Keck spectra plus the single template spectrum, computed as de- 
scribed in the text. The mean value has been subtracted. Vertical 
scale for all three panels is the same. 

constraint on eccentricity. We found an eccentricity of 
e = 0.046 ± 0.031, not significantly different from zero. 
Therefore the data wore rc-fit with a circular orbit, in 
which both the period P and the mid-transit time Tc 
were held fixed at their values from the light curve anal- 
ysis (Table 3). The solution fits the data well; see Fig. 2a. 
No long-term trends are seen in the residuals. 

In order to get a reduced chi-squarc value near unity 
for the fit, it was necessary to add an additional random 
noise component with amplitude 8.6 m in quadrature 
to the internal errors. This is essentially identical to the 
8.7 ms""'^ radial velocity "jitter" expected to arise from 
stellar surface activity, based on the above-mentioned 
strength of the emission cores of the Ca+ H and K lines, 
^ogR'HK = -4-81 (Wright 2005). The parameters of 
the resulting final fit are not significantly changed by the 
inclusion of this jitter, and are listed in Table 3. 

Following Torres ct al. (2007), we explored the pos- 
sibility that the measured radial velocities arc not real, 
but instead due to distortions in the spectral line profiles 
due to contamination from a nearby unresolved eclipsing 
binary. In that case the "bisector span" of the aver- 
age spectral line should vary periodically with amplitude 
and phase similar to the measured velocities themselves 
(Queloz et al. 2001; Mandushev et al. 2005). Instead, we 
detect no variation in excess of the measurement uncer- 
tainties (see Fig. 2, bottom panel). We conclude that the 
velocity variations are real and that the star is orbited 
by a Jovian planet. 

Combined with the results from the light curve and 
radial velocity modeling, the above stellar parameters 
yield a planet mass of Mp = 1.06 ± 0.12 Afj. The surface 
gravity, gp = 14.8 ± 1.0 ms^"'^, was obtained from the 
light curve and radial velocity fits (see Southworth et 
al. 2007). The mean density and its uncertainty, Pp = 
0.558 ± 0.047 gcm~^, follow from the absolute radius of 
the planet. All planet parameters are listed in Table 3. 

5. DISCUSSION 

HAT-P-6b, with radius 1.33 is similar in size to 
five low density "inflated" planets tabulated by Kovacs 



ct al. (2007) (i.e. WASP-lb, HAT-P-4b, HD 209458b, 
TrES-4, and HAT-P-lb). However, its mass of 1.06 Afj 
is greater than the mass of any of these, and hence it has 
a larger mean density and surface gravity. For a planet 
of its mass, age of 2.3 Gy, and stellar flux Fp at the 
planet given by Fp = L^/ {Ana^), models of Burrows et al. 
(2007) predict a radius of about 1.21 Rj, assuming that 
the planet has no heavy-element core. The metallicity of 
HAT-P-6, [Fe/H] = -0.13±0.08, is among the smallest of 
known transiting planet host stars. If we assume that the 
bulk composition of IIAT-P-6b tracks the metallicity of 
its host star, the size of its heavy-element core should be 
small, but not vanishingly so. Thus the predicted radius 
from Burrows et al. (2007) is comparable to, but perhaps 
slightly higher, than one might expect for IIAT-P-6b. 
However, the actual radius found here, Rp = 1.330 ± 
0.061, lies above the predicted value by about 2a, so it 
appears to be somewhat inflated relative to that model. 

Hansen & Barman (2007) proposed that hot jupiters 
can be placed into two classes based on their equilib- 
rium temperature and Safronov number 0, where @ = 
(a/Rp) X {Mp/Mi,) is the ratio of the escape velocity 
from the surface of the planet to the orbital velocity. 
When Safronov number is plotted versus equilibrium 
temperature, transiting hot jupiters seem to fall into two 
groups, with an absence of objects between. However, 
the Safronov number for HAT-P-6b is 0.064±0.004; along 
with HAT-P-5b (Bakos et al. 2007b) with Safronov num- 
ber 0.059 ± 0.005, these two planets appear to fall be- 
tween the two groups in such a plot. Hansen & Barman 
also noted a difference in the relation between planet 
mass and equilibrium temperature for planets of the two 
classes, but HAT-P-6b and HAT-P-5b appear to fall be- 
tween the two classes in this respect as well. It would 
seem that discovery and characterization of a large num- 
ber of additional transiting exoplanets may be neces- 
sary to establish unambiguously whether there is a bi- 
modal distribution of hot jupiter planets according to 
their Safronov number. 
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